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Abstract 

We obtain analytic formulae for the null geodesies of Friedmann-Lemaitre- Robertson- 
Walker spacetimes with scalar perturbations in the longitudinal gauge. We use these to 
provide a rigorous derivation of the cosmological lens equation. We obtain an expression for 
the magnification of a bundle of light rays in these spacetimes without restriction to static 
or thin lens scenarios. We show how the usual magnification matrix naturally emerges in 
the appropriate limits. 
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ly-j ■ 1. Introduction 

o . 

The bending of light by a single symmetric gravitational lens in a Euclidean space is 
shown in Fig. 1. The symmetry guarantees that the lines of sight from the observer, o, 
to the lens, I, and to both the lensed, and unlensed, image of the emitter, e, lie in the 
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same plane, as do the angles a, f3, and 9. D(o,l) and D(o,e) are the distances from the 
observer to the lens and emitter, respectively. We assume that the deflection angle, a, is 
small. Then, locally about the line of sight to the emitter's image, we may approximate 
the two-spheres at distances D(o,l) and D(o,e) from the observer as planes, called the 
lens plane and source plane respectively. D(l, e) denotes the distance between these two 
^ ■ planes. The assumption of small deflection angle also allows us to relate the lensing angles 
by (3 = 9—aD(l, e)/D(o, e). This is the simplest example of the gravitational lens equation. 

The generalization of this equation to more complicated lens structures and non- 
Euclidean background spaces proceeds by a number of steps. A general lens is not sym- 
metric so that the angles a, f3, and 9 are not necessarily coplanar. To handle this, consider 
a set of cartesian axes with origin at the observer. Choose the a>axis to coincide with 
the line-of-sight to the image. Let a 1 , where % runs over {2, 3}, be the angle between the 
x-axis and the projection into the x l x-plane of the line-of-sight vector from the deflection 
point, p, to the emitter. Similarly, let j3 l be the angle between the projection into the 
x l x-plane of the line-of-sight vector from the observer to the lens and the projection into 
the x l x-plane of the line-of-sight vector from the observer to the unlensed image. Also, 
let 9 l be the angle between the x-axis and the projection into the a^x-plane of the line-of- 
sight vector from the observer to the lens. To allow for non-Euclidean spatial geometries 
D(o, e) is taken to be the angular-diameter distance in the background geometry from the 
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observer to the intersection of the x-axis with the source plane, and D(l, e) is taken to be 
the angular-diameter distance in the background geometry between the deflection point, 
p, and the intersection of the x-axis with the source plane, p' . Then, again assuming small 
deflection angle, 



pi = ei _ (i) 

This is the standard cosmological gravitational lens equation (Schneider, Ehlers, and Falco 
1993). The magnification matrix, M l j = d(3 l /dd^ contains information about the defor- 
mation of ray bundles connecting the observer and emitter. For example, the inverse of 
the determinant of this matrix is the magnification of an image relative to an unlensed 
image. 

The purpose of the current paper is to address a number of subtle issues that arise 
when we attempt to justify mathematically the use of the lens equation (1) for calcula- 
tions in our Universe, although most workers agree that the physical justifications for its 
use in observed lens systems are strong. First, there is the question of the best choice of 
distance factors. There exists a large literature addressing this question, primarily con- 
cerned with the appropriateness of the so-called Dyer-Roeder distances (Dyer and Roeder 
1972, 1973; Ehlers and Schneider 1986; Futamase and Sasaki 1989; Watanabe and Tomita 
1990; Watanabe, Sasaki, and Tomita 1992; Sasaki 1993). Our work suggests that within 
the framework of cosmological perturbation theory, the natural distance factors to use are 
those of the background. Hence, the choice of distance factors is equivalent to the choice 
of cosmological model, in agreement with the recent results of Sasaki (1993). The issue of 
which cosmological model is most appropriate must be addressed in its own right. 

The second issue which should be addressed in any mathematical investigation into 
the lens equation concerns the accuracy of the approximation of an actual photon path by 
two geodesies of the background which join at a point near the lens, the deflection point, 
p. On physical grounds we expect this approximation to be good for systems for which 
the photon-lens interaction is localized: the thin-lens approximation. One purpose of our 
present work is to quantify the relationship between the actual path and that used in the 
lens equation. 

Third, how are we to find the angles appearing in the lens equation from physical 
data? Generally, the a 1 are taken to be those predicted by calculations in Einstein-de 
Sitter spacetime, since the overall curvature of space should not be important near p, where 
the light ray interacts with the lensing object. For static, thin lenses, these calculations 
write the deflection angle as a superposition of point mass deflection angles contributed by 
each mass element of the lens projected onto the lens plane (Schneider et al. 1993). For 
brevity we will term the resultant angle the "Einstein angle." Another purpose of this work 
is to derive this result from the full equations of light propagation under an appropriate 
set of mathematical approximations. The lens equation, for static, thin lenses, effectively 
assumes that the light path is described by the Jacobi equation of the background spacetime 
subject to an impulsive wavevector deflection at the lens plane by an angle equal to the 
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usual Einstein bending angle. We wish to quantify the level of approximation involved. 

There have been two notable recent attempts to clarify the validity of the cosmological 
lens equation by deriving it from the optical scalar equations (Seitz, Schneider, and Ehlers 
1994) and the Jacobi equation (Sasaki 1993). However, a crucial difference between these 
papers and the present work is that they treat the path of the light ray differently near to 
and far from the lens. It is precisely this assumption that we must eliminate if we hope to 
gain a more general lens equation able to quantify the errors implicit in equation (1). 

Our approach is to investigate the cosmological lens equation as it emerges from the 
geodesic equation. We are able to do this by making use of a technique for constructing 
null geodesies in perturbed spacetimes (Pyne and Birkinshaw 1993). The conventions of 
the present paper are the same as those in that earlier work. The results of this letter come 
from applying this technique to FRW spacetimes with scalar perturbations in the longitu- 
dinal gauge. The calculations yielding the results presented here will appear in Pyne and 
Birkinshaw (1994). The results we report here are: analytic formulae (equations (4) and 
(5)) for light rays in the spacetime (2); a general expression for the magnification under- 
gone by a bundle of light rays capable of handling non-static, geometrically thick, lenses 
(equations (10), (11), and (13)); and a demonstration that the usual deflection angle, lens 
equation (1) and magnification matrix are recovered in the appropriate approximations. 
To our knowledge, these are the first rigorous demonstrations of these results for perturbed 
FRW spacetimes. 

2. The Deflection Angle 

Our starting point is a choice for the metric of our Universe. We choose to work with 
FRW spacetimes with scalar perturbations written in the longitudinal gauge, 

ds 2 = a 2 [-(1 + 2(P)dr] 2 + (1 - 2(P)~f- 2 (dx 2 + dy 2 + dz 2 )] (2) 

where 7 = l+Kr 2 /4, k being the spatial curvature parameter (±1 or 0) and r 2 = x 2 +y 2 +z 2 . 
Inhomogeneities are represented by the quasi- Newtonian potential, 0. For the order needed 
by us the expansion factor, a, is unperturbed from its Friedmann form (Jacobs, Linder, 
and Wagoner 1993). The metric (2) is also used by Seitz et al. (1994) and Sasaki (1993), 
so that direct comparison of results is possible, and recent work by Futamase (1989) and 
Jacobs et al. (1993) has shown that structure of galactic scales and greater in our Universe 
can be well modeled by metrics of this type. While results obtained with the metric (2) are 
not appropriate for lensing by gravitational waves or vector perturbations, it is possible to 
treat these cases by a similar method (see Pyne and Birkinshaw 1994). 

Because the Friedmann expansion, a, plays the role of a conformal factor it is simplest 
to work with the null geodesies of ds 2 , defined by ds 2 = a 2 ds 2 . Quantities in ds 2 will 
always be written with an overbar. Light rays in these two metrics coincide and their 
(affine) parameterizations are related by k^ = a~ 2 k^. While our formulae are mostly 
written in terms of quantities in ds 2 , our results apply to the actual spacetime, ds 2 . We 
note that by locating the observer at the spatial origin the radial null geodesies of ds^ 2 , 
(i.e. of that part of ds 2 independent of <f>), obey k^° = 1 and k^ 1 = —7c 1 , where e 1 are 
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the direction cosines at the observer, so that Xli=i( e *) 2 = 1 (McVittie 1964). The explicit 
solutions for the comoving radius and for 7 along such rays are given by 



r(A) = 2 tan K ' 



7(A) = sec K 



2 / A — A 



(3) 



where A D is the affine parameter at the observer. The subscript k on a trigonometric 
function denotes a set of three functions: for k = 1 the trigonometric function itself, for 
k = — 1 the corresponding hyperbolic function, and for k = the first term in the series 
expansion of the function. The paths of the rays are = A, x^ 1 = re 1 . 

The method presented in Pyne and Birkinshaw (1993) allows us to write the null 
geodesies of ds 2 as x^ = x^ ^ + x^'^. For light rays passing through the observer 
given as above and the separation vector is 



x y 



!)°(A) = +2 / du(u - \)k^ m (j) >m {u) 



(4) 



and 



x 



(i)< 



(A) = -2k Wl [ du(u-X)^(u) 
Jx drj 



+ 2 7 (A) / dusm K (u- A)^-(V ± (/)) , 



(5) 



where V ± % = ^(°) mn (^ — e n e l ) d m is the transverse gradient operator. We see that the 
spatial separation is naturally written as a the sum of a longitudinal and a transverse 
term. The above integrals are taken over the background geodesic Hence if these 

solutions are used for geometrically thick lenses, it will in general be necessary to apply 
an iterative procedure, incorporating a number of background paths. Such an approach is 
familiar from the usual multiple-lens plane theory (Schneider et. al. 1993) but we stress 
that formulae (4) and (5) allow the photon path be approximated to arbitrary accuracy, in 
contrast with the multiple lens plane theory where the continuum limit is not compatible 
with the asumptions underlying the theory. Precise statements of consistency will appear 
in Pyne and Birkinshaw (1994). 

We can understand (4) and (5) by considering their relation to the equation of geodesic 
deviation. The Jacobi equation of ds^ 2 subject to an arbitrary impulsive wavevector 
perturbation dk^ at some affine parameter u, is solved by deviation vector 5x^ with spatial 
components 
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5x l (\) = sm K (u - \)5k l ± (u) - lQ(u-\)6kUu) (6) 

where 5k]_ = (5] — e l ej) bW is the impulse in the transverse direction, and 8k^ = 5k l — 8k\ 
is the longitudinal impulse. 

Comparing (5) with (6), produces the following interpretation of the spatial compo- 
nents of our solution for the separation, x^ % : the photon path differs from a background 
path because of a continuous sequence of impulsive perturbations 

Sk* = -2 (V ± <PY + 2£; (0) ^. (7) 

In fact, the form of the impulse can be gained directly from our equations by differentiating 
the spatial separations with respect to the affine parameter, and inserting a delta function 
at A; into the integrand which forces the integrand to vanish except at the lens plane. This 
gives 

fcW (AO = -2 (V±0)' (Ai) + 2fc<°>' (A/) ^ (A;) , (8) 

which is exactly the impulse found above, equation (7). 

At this point we need only do a little work to recover the Einstein angle from our 
equations. Establish a set of cartesian axes at the observer and choose the unperturbed 
wavevector k^ = (1,-7,0,0). Consider a static, localized perturbation in the xy-plane. 
Then the angle represented by the impulse perturbation found above, that is, the angle 
k^ 1 (A/) + k^ 1 (A/) makes with k^ 1 (A/), is given by 

* V = \fr =-2 7 (W,„. (9) 

The factor of 7 (A;) is present only because our co-ordinates are scaled in an unusual way 
at the lens plane. If we make the co-ordinate change x 1 ' = x t / , y(Xi) the metric on the 
lens plane becomes Minkowskian and, locally near the deflector, y' and z' serve as normal 
co-ordinates on the lens plane. In these co-ordinates on the lens plane takes on the 
usual Newtonian form (with origin shifted from the lens, accounting for the unusual minus 
sign). Since our lensing angle — 27(A/)0 j2/ = —20^/ the integrated impulse lensing angle 
for a localized perturbation is exactly the Einstein deflection angle. We emphasize that 
this is the first time this has been shown rigorously for the curved FRW spacetimes. 

For completeness we note that the timelike component of the separation can also 
be analyzed by comparison to the Jacobi equation. The Jacobi equation of ds^ 2 for an 
impulse wavevector perturbation 5k^ at affine parameter u results in a timelike component 
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of the deviation vector 5x° = — (u — X)Sk°(u). Comparison with our solution for the 
separation reveals that the time delay undergone by the light ray relative to the fiducial 
background ray may be considered to result from a sequence of impulses 5k° = —2k^ m (f) jTn . 

3. The Magnification 

We now want to examine the magnification undergone by a bundle of light rays. We 
define this after Schneider et al. (1993) in the following way. Suppose a source of given 
physical size at some redshift is observed to subtend solid angle dfl. An identical source 
observed at identical redshift placed in an FRW spacetime would subtend solid angle dQ(°\ 
The magnification M is defined to be dfl/dfl^°K 

We can construct a bundle of light rays in the spacetime (2) which emanate from a 
source and converge at an observer by varying the direction cosines, e l , in equations (4) 
and (5). This yields dQ in terms of the redshift and proper size of the source. A similar 
relationship is easy to derive for a source in the background spacetime. Comparison of 
the two yields M. If the four-velocities of the observer and emitter are written u^ e ^ = 

(1 — o ( e ), v l t e )), which are correct to first order in (j) and v, then 



M = (l + 20 o + 2 
/ (DetMy 



v*ki 0) 



- 2fcW° + 2 cot, (A - A e ) SX e + k sin, (A c - A e ) e t x^ 1 ) 

o / 



(10) 

where [f] e Q = f (A e ) — / (A D ) and the subscript e indicates evaluation at A e . Here 5X e 
enforces the equal redshift aspect of the definition. It is set by l + z(X e +5X e ) = l + z(°)(A e ) 
and is given by 



6X e = ^ 
a e 



^(0)_aa- 1 xW° + ^ 1 )°] e 



(11) 



M l j is the magnification matrix, of dimension 2x2, 



M\ =I 2 + - 



1 dx^ 1 



r dei 



(12) 



where I2 is the two-dimensional identity matrix and the indices i,j run only over the 
transverse spatial dimensions. For simplicity, we will henceforth assume the unperturbed 
wavevector fc(°) = (-1, 7 ,0,0) so that ij run over {2, 3}. Keeping in mind the transverse 
nature of these indices we can write d/de^ — djdd\ with 9 the vectorial angle of equation 

(!)• 

Explicit calculation yields 
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2 f Xe 36 

3 sm K (A - A e ) 3 J Xo dr] 

2 f Xe 

- 7T 7-r / rf« sin K (u - A e ) 7(it) (13) 

sin* (A - A e ) J Xo 

x Si^-riuh-^g^iu^iu) 

To see how the usual magnification matrix is contained in equation (13) we rewrite the 
second integral using the angular diameter distance of ds^°\ D(u, A) = a(A) sin K (w — A). 
In this way 



with the two-dimensional projected angle 



da 1 



(14) 



a = -2 — -— . (15) 

7(u) 

We now consider lensing by a static, geometrically thin lens. The first term in equation 
(14) vanishes and we approximate the angular diameter factor as constant over the region 
for which the potential is important so that 



D(Xi,e) d i; 
D(o, e) dOi J Xo 



^«*S---n^^7/ dua* (16) 



We have already seen that the integral of 6? over the background path produces the 
Einstein deflection angle. As a result we conclude that our equation has reproduced the 
usual magnification matrix defined by the 9- gradient of equation (1). 

4. Summary 

We have presented formulae (4), (5) for the null geodesies intersecting an observer's 
worldline in an important class of perturbed spacetimes, FRW backgrounds with scalar 
perturbations, in the longitudinal gauge. We have used these equations to obtain a general 
formula (10) for the magnification of ray bundles in these spacetimes. With this, we can 
show for the first time how the usual lens equation (1) and magnification matrix are 
recovered in these spacetimes without dividing light paths into near and far lens regions. 
In forthcoming papers we will consider the implications of these results for practical lensing 
calculations. 
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Figure Captions 

Figure 1. Gravitational lensing by a single symmetric lens. 
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